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Abstract 

We consider the ground state of simple quantum systems coupled to an environment. In general the system is entangled with 
its environment. As a consequence, even at zero temperature, the energy of the system is not sharp: a projective measurement 
can find the system in an excited state. We show that energy fluctuation measurements at zero temperature provide 
entanglement information. For two-state systems which exhibit a persistent current in the ground state, energy fluctuations 
and persistent current fluctuations are closely related. The harmonic oscillator serves to illustrate energy fluctuations in a 
system with an infinite number of states. In addition to the energy distribution we discuss the energy-energy time-correlation 
function in the zero-temperature limit. 
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1. Introduction 

A quantum system cooled to zero temperature nev- 
ertheless knows about its environment since generically 
the system state and the bath state are entangled. The 
ground state does not factorize into a product of a sys- 
tem wave function and a bath wave function. Entangle- 
ment of two subsystems [1] is often discussed in terms 
of the strange non-local properties it implies for sys- 
tems that can be spatially separated. Here we consider 
two-level systems, often now called qubits, or harmonic 
oscillators which are coupled to a bath. In such a "ther- 
modynamic" setting we can not easily separate the two 
systems and apply a Bell test [2] to verify the entan- 
glement. Nevertheless, systems entangled with reser- 
voir states exhibit a number of properties which distin- 



^ Corresponding author. 

E-mail adress: Markus.Buttiker@physics.umge.ch 



guishes them from systems for which the ground state 
factorizes. 

A particularly instructive quantity to consider is the 
energy of the system. First, from a purely theoretical 
point of view, we always have the energy of the sys- 
tem as an observable. We write the energy of the total 
system as 

H = Hs + Hc + Hb (1) 

where Hs is the system. He the coupling and Hb the 
bath energy operator. A second, more important rea- 
son for considering the energy is that at zero tempera- 
ture fluctuations in energy are a direct indicator for sys- 
tem bath entanglement [3]. In contrast, if the system- 
bath state is not entangled, the system is simply in its 
lowest energy state. This case is often viewed as self- 
evident instead of the generic case addressed here. 

Consider a two level system with energies E+ and E- 
and probabilities p± to find the system in the excited 
level and in its ground state. The expectation value of 
the energy is 
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Coupling constant a 

Fig. 1. Energy probabilities p+ and p_ for the spin-boson 
problem as a function of coupling strength a and different e. 
For weak coupling to the environment increasing coupling 
increases the the probability to find the qubit in the excited 
state, with increasing coupling to the environment it it is 
more likely to measure the qubit in the excited state. After 
Ref.[3] 
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and the energy fluctuation away from its average 

For a moment consider the system to be in ther- 
mal equilibrium coupled to a bath at temperature 
T. Later we will only consider the zero-temperature 
limit. Assume for the moment that there is no en- 
tanglement. Then statistical mechanics tells us that 
p+ oc exp(— i5+/fcr) and p_ oc exp{—E-/kT) or with 
proper normalization (p+ +p- ~ 1) 

1 



P± 



(4) 



l+exp{±(E+ - E^)/kT)' 

As the temperature tends to zero, p_ tends to one (p-i- 
tends to zero) and consequently the energy fluctuations 
Eq. (3) tend to zero. The energy of the system is a 
sharp variable in the ground state. 

This simple picture is dramatically modified if the 
system and the bath are entangled. Now statistical me- 
chanics must be applied to the entire system (system 
plus bath plus interaction energy). The system alone 
is now determined by a (reduced) density matrix with 
elements p±,± on the diagonal and non-diagonal ele- 
ments p±,^. In the energy eigen basis of the two state 
system, p± = p±,±- While the ground state of the en- 
tire system is in general a pure state, the density matrix 



of the system alone is in general that of a mixed state. 
For the simple model of a two state system (a spin) 
coupled to a harmonic oscillator bath (the spin-boson 
problem [4,5]) an ohmic bath leads for weak coupling 
to a probability [3] 



p+ = Qlog(a;c/n). 



(5) 



Here a is the system bath coupling constant, hQ — 
E+ — E- is the energy separation of the two levels and 
u!c is a cut-off of the oscillator spectrum of the bath 
(a Debye frequency). As a consequence, even at zero 
temperature the energy of the system is not sharp but 
fluctuates according to Eq. (3) 



(6) 



We emphasize that this result is not special for a 
two state system. For instance for a harmonic oscil- 
lator with frequency uio coupled to an ohmic bath 
of harmonic oscillators, we find that the energy 
fiuctuations[6,7] are given by 



(7) 



Thus huo plays a role similar to the energy separation 
hQ of the two level system. 

The energy distribution is certainly far from being 
Gaussian. Thus the mean square deviations given by 
Eqs. (6) and (7) might not be a good indicator of the 
way the energy is distributed over the different states 
of the system. For instance for the two state system 
discussed above the probability distribution P{E) to 
find find the system with energy E in the interval dE 
obviously consists of two deita-function peaks at E+ 
and E~. Thus we can write this distribution in the 
form, 

P{E)^p+S{E-E+)+p-5{E-E-) (8) 

where p± are as above the probabilities to find the 
system in the excited state and in the lowest energy 
state of the system. Now since p_ 4- p+ = 1 and since 
{E) = (Hs) = E-p- + E+P+ we can also express the 
probabilities p± in the form p+ — {E)/hQ — E-/hfl 
and p- — E+/hfl — {E)/hQ. Here we have used hQ — 
E+ — E- . Without loss of generality we can set E± — 
±hQ/2 and thus flnd for the distribution [3] 

(9) 
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Note that if the system and bath are decoupled we have 
(E) — —hQ/2 and the distribution function consists of 
only one peak at i5 = —hQ/2. 

We emphasize that (E) = (Hs) is the expectation 
value of the system's energy in the overall ground state 
of the system plus bath plus interaction energy. Below 
we will also discuss P{E) for the harmonic oscillator. 
For a harmonic oscillator P{E) consists of an infinite 
number of deita-functions with a rapidly decreasing 
weight of the higher lying states. 

To be specific, consider the two-state system Hamil- 
tonian to be 



kT* = 



H, = (e/2) a, + (A/2) 



(10) 



where e measures the distance from resonance and A 
is the energy separation at resonance, az and are 
Pauli spin matrices. The level separation is thus de- 
termined by the frequency fl = \/e^ + A'^/h. If this 
system is coupled via to a harmonic oscillator bath 
(spin-boson problem) Hs will not commute with the 
total Hamiltonian. For this system the probabilities to 
find the system in the excited state and in the low en- 
ergy state are shown in Fig. (1). At resonance e = 0, 
the probabilities p± tend with increasing coupling con- 
stant towards 1/2. For these parameters we can use the 
Bethe solution of the anisotropic Kondo model [8,9]. In 
the strong coupling limit an energy measurement will 
find the system with equal probability in the ground 
state and in the excited state of the system. If the sys- 
tem is not symmetric e 7^ the probability to find the 
system in an excited state reaches a maximum as a 
function of a and tends to zero for very strong cou- 
pling. For large e the probabilities can be found from 
perturbation theory [8,9]. Combining both the Bethe 
Ansatz solutions and the perturbation theory [9] it is 
possible to give the probabilities over the entire range 
of parameters. Computational work based on renor- 
malization is reported in Ref . [10] . 

Experiments are always carried out at finite temper- 
ature, and it is important to demonstrate that there 
exists a cross-over temperature to the quantum behav- 
ior discussed here. In the low temperature limit, the 



thermal occupation probability is p+ = e' 



-nn/kT 



(see 



Eq. (4). In the weak coupling limit for the symmet- 
ric spin boson problem, the probability to measure the 
state as "spin up" scales as p+ = a\og{LOc/^)- Setting 
these factors equal and solving for T* yields [3] 



log(alog^) 



(11) 



Since T* scales as the inverse logarithm of the cou- 
pling constant, it is experimentally possible to reach a 
regime whore thermal excitation are negligible. Exper- 
imentally and theoretically [11] one might be tempted 
to define temperature with the help of the qubit by 
fitting p+ = e^'^^^^^'^''ff with an effective temperature 
to the experimental data. Experimentally one would 
then find that it is impossible to cool the qubit below 
the temperature T*. Still a careful examination would 
show that the state is not in fact "thermal" since it 
depends on the coupling constant to the bath. 

The temperature T* can be viewed as a measure of 
the energy difference between the energy of the lowest 
energy separable state and the true entangled ground 
state. Recent works emphasize the role of this energy 
difference as an entanglement witness [12,13]. We now 
relate energy fluctuations directly to entanglement in 
the ground state. 



2. Energy fluctuations as an entanglement 
witness 

Energy fluctuations arc determined by probabilities 
alone, that is by the diagonal matrix elements of the 
density matrix only. Thus it is not obvious that we 
can make a statement about entanglement. In general 
such a statement also depends on the non-diagonal el- 
ements of the density matrix. However, we are given 
the additional information that we are in the ground 
state. If we measure the subsystem's energy and find 
an excited energy, then we know that the state is en- 
tangled. For weak coupling we can make quantitative 
statement. The probability to find the system in higher 
lying states is exponentially suppressed. To first order 
in the coupling constant, we can consider a two-state 
system where the density matrix has the form p- = 
p — = 1 - op, p+ = P++ = ap, p+- = p*_+ = ac. 













p = 




+ a 
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For vanishing coupling constant q = 0, this just gives 
the density matrix for the separable state. The linear 
dependence of p on a holds to first order for the model 
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systems considered here and is the entanglement con- 
tribution. If one measures the diagonal elements of p, 
one obtains p — and p++ as the probability to measure 
the system in the ground or excited state. If we now di- 
agonalize p, the eigenvalues are \± = {1 — p a,pa} -\- 
0{a^). To first order in Q, the eigenvalues are the di- 
agonal matrix elements, so we may (to a good approx- 
imation) write entanglement measures like the purity 
in terms of these probabilities. As an example the pu- 
rity Tr(p^) is given by Tr(p^) = A+ + Ai ~ 1 - 2ap = 
l-2p+. 

We next consider two examples of qubits which ex- 
hibit the behavior discussed above. 



3. Persistent current qubits 

Thus far we have focused on the energy of the system 
as the quantity of interest. In this section we show that 
other observables which reflect properties of the system 
behave in fact very similarly. In particular we discuss 
two states systems (qubits) which have the property 
that their tunnel matrix element A in Eq. (10) is de- 
pendent on an Aharonov-Bohm flux 3>. The free energy 
of such a system then depends on the flux and in gen- 
oral supports a persistent current in its ground state 
/($) = —dF/d^. The prediction of persistent currents 
in small normal and disordered loops [14] has played 
an important role in the development of mesoscopic 
physics and continuous to be a subject of current in- 
terest [16]. 

3.1. The mesoscopic persistent current qubit 

A small metallic loop shown in Fig. (2) can be made 
into a two state system with the help of a quantum 
dot [17]. For sufficiently small charging energy only the 
states with N and N + 1 electrons on the dot will be 
relevant. As long as these two states are energetically 
very different, charge on the dot is fixed and transport 
through the dot is blocked. Only near the point of de- 
generacy can electrons tunnel in and out of the dot. 
In the charge basis, the charge on the dot is propor- 
tional to (Tz with an energy e that determines how far 
away the system is from the point of degeneracy. At 
the point of degeneracy, charge tunneling is permitted 
since the dot is weakly coupled with tunnel energies tL 




Fig. 2. Ring with an in-line quantum dot penetrated by an 
Aharonov-Bohm flux. The ring is coupled capacitively to 
an external circuit. After Ref. [20] 

and tn to its contacts. In such a ring the persistent cur- 
rent exhibits sharp peaks at special values of the gate 
voltage [17] much like the Coulomb peaks of conduc- 
tance. The Aharonov-Bohm flux can be incorporated 
into an effective tunnel matrix element which connects 
the charge states A'^ and N + I, 

AV4 ^tl + t%± 2tLtR cos(27r<I>/$o). (13) 

Here the sign depends on the number of electrons in the 
ring. Thus near a point of degeneracy the Hamiltonian 
of this ring is of the form given by Eq. (10). We refer 
the reader to Refs. [17] for a detailed derivation. 

The system is coupled to an external circuit via ca- 
pacitances. In particular, if the external circuit is an 
ohmic resistor, we can replace it with a transmission 
line with equal input impedance [9]. The transmission 
line represents a harmonic oscillator bath [18] and the 
entire system is a particular realization of the spin- 
boson problem with the interesting feature that the 
tunnel matrix-element is flux dependent. 

The model permits us to address the interesting 
question of how persistent currents are affected by 
environments. We follow here the discussion of Ce- 
draschi, Ponomarenko and one of the authors [9]. An 
extended discussion of the initial work is provided in 
Ref. [19]. The persistent current is obtained using the 
results of Bethe Ansatz and perturbation results for 
the anisotropic Kondo model [8]. For weak system- 
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bath coupling, Ref. [20] provides a discussion based on 

a quantum Langevin approach. Closely related works 
investigate the effect of a fluctuating Aharonov-Bohm 
flux [21], the effect of charge fluctuations in nearby 
pure and disordered conductors [22,23] and the effect 
of hot bosonic baths [24]. 

Let us now first show that the persistent current, 
like the energy of the system, is not sharp but fluctu- 
ates. We are interested in current fluctuations so we 
need an expression for the current operator. Since the 
persistent current is due to electrons of the ring alone 
we can consider the isolated system a = 0. Eq.(lO) is 
the Hamiltonian in the charge basis. In the eigen ba- 
sis the Hamiltonian is simply Hs = (li.fl/2)az. In the 
eigen basis the persistent current carried by a state 
is determined by the derivative of the energy of this 
state with respect to flux. The persistent current of 
the eigenstates with energies ±hfl/2 is =fJo with Jo = 
{l/2)dMl/d^. One of the states corresponds to a clock- 
wise persistent current and one state corresponds to a 
counter-clockwise current. Thus in the eigen basis the 
current operator is simply 

/ = /ocT^ = -{2h/n^)Hs. (14) 

Here we use that in the eigen basis, both the current 

operator and the energy arc proportional to o^- Thus 
the persistent current is directly related to the energy 
of the system [25]. While the first expression in Eq. 
(14) is valid only in the energy eigen state basis, the 
second expression is in fact general. 

In the presence of a bath, the persistent current, like 
the energy of the two state system, is in fact not sharp 
but fiuctuates. Since the average persistent current is 
(7) = I+P+ + I-p- with p+ -h p- = 1 we can also 
express p± in terms of the persistent current. With 
7± = ±/o we find p± = (1/2)(1 =f {i)/h) and the 
distribution of currents is thus 

P{I) = i(l + T^)5(/-/o) + i(l-^)<5(/+/o) (15) 

Note that p± = (1/2)(1 =F (I)/Io) < 1 implies that 
coupling the ring to the bath can only suppress [9] the 
persistent current {/) < /q. 

Next let us investigate the flux dependence in more 
detail. In the absence of coupling to a bath a = 0, and 
at the point of degeneracy e = the persistent current 
is, 7o($) = dA($)/d3>. Consider now additionally the 
special symmetric case, when the tunneling rates are 




Fig. 3. Persistent current as a function of the 
AhaJTonov-Bohm flux at resonance (e = 0) for different cou- 
pling strengths a. Arbitrary small coupling suppresses the 
discontinuity of the a = 0-persistent current at $/$o = 
and 1. 

equal t = tL = ifl. Taking the lower sign in Eq. (13) the 
persistent current is 7($) = — (e//i)47ri cos(7r$/<l>o) in 
the interval < 3> < 3>o ■ It is a periodic function with 
a discontimious jump at <I> = n$o, n = 0, ±1, ±2, ... It 
is shown as a solid line in Fig. (3). 

It is now very interesting to investigate what hap- 
pens to the discontinuous jump in the persistent cur- 
rent in the presence of the bath. Since a Fourier repre- 
sentation of the persistent current 

7($) = ^ 7„ sin(27rn$/*o) (16) 

n 

needs arbitrary high harmonics such a jump is a sig- 
nature of a perfectly coherent system: the top-most 
electron in our system which gives rise to this current 
must circulate the Aharonov-Bohm flux coherently n- 
times to generate the n — th harmonic. Coupling to a 
bath generates de-coherence, and we suspect that the 
bath suppresses such a discontinuity immediately. The 
Bethe Ansatz solution of the spin-boson model gives a 
persistent current 

/($) tx A'T^^dA/d$ (17) 

where a is the coupling constant. The exact result is 
given by Cedraschi et al. [9]. The persistent current 
as a function of flux for different coupling constants is 
shown in Fig. (3) . Thus an arbitrary small coupling to 
a bath is sufficient to suppress the discontinuity in the 
persistent current. Using Eq. (17) a more quantitative 
analysis can be provided. Pilgram [26] finds for the 
Fourier amplitudes 
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Fig. 4. Cooper pair box: A superconducting metallic dot 
is coupled via two Josephson junctions (with energy Ej 
and capacitance Cj) to a superconductor terminal and is 
capacitively coupled (with capacitance Cg) to a gate per- 
mitting the control of charge. An external flux through 
the hole of the structure controls the Josephson energy. 

J„ _ (2Q-l)...((2n-2)a-(2n-3) 

7i (4a - 5) (2na - (2n + 1)) ' 

With the Ansatz 

^=An exp(-6„(Q)a(n - 1)) (19) 

the values found for 6„(a) are, &2 = 6/5, 63 = 88/105, 
64 = 626/945... Thus the bath suppresses the persis- 
tent current almost in an exponential manner, as if the 
system were subject to dephasing. This suppression is 
the stronger the higher the harmonics. 

The system of Fig. (2) can viewed as a double quan- 
tum dot. Recent experimental work by Hayashi et al. 
[27] has demonstrated that conditions can be achieved 
for which the double dot system is well described by 
a Hamiltonian of the form Eq. (10). The work also 
demonstrates how strong measurements can be imple- 
mented. Refs. [28,29] represent closely related theoret- 
ical work. 

3.2. The split Cooper pair box qubit 

In this section wc compare briefly a superconduct- 
ing structure [30,31,32,33] with a behavior that is anal- 
ogous to the model discussed above. This example is 



important since, in contrast to the normal state qubit 

discussed above, it does not depend on single parti- 
cle energies and since a projective measurement test- 
ing the state of the system by measuring its persistent 
current, has in fact been implemented [33]. A review of 
the current research on superconducting qubits is pro- 
vided by Devoret, WalrafFand Martinis [34]. The struc- 
ture of interest hero is a Cooper pair box [35] , a small 
superconducting island coupled with tunnel junctions 
to a large superconductor and coupled capacitively to 
a gate (see Fig. 4). The small superconducting island 
can be split [32] such that it forms together with the 
large superconductor a ring. If the capacitances of the 
junctions Cj and the capacitance to the gate Cg are 
small, the energy for charging the island with an ad- 
ditional Copper pair Ec — (2e)^/(2Cj -I- Cg) is large 
and the system can effectively bo described with a two- 
state Hamiltonian. The limit of interest is the charge 
controlled Cooper pair box [34] for which the Joseph- 
son energy Ej is much smaller than the charging en- 
ergy Ec- The parameters [34] of the two-state Hamil- 
tonian Eq. (10) are e — Ej cos{tt^x/^o) where $2, is 
the externally applied flux, <I>o = h/2e the charge 2e- 
flux quantum and A = Ec ( 1 /2 — Ng ) where Ng can be 
controlled by adjusting a gate voltage. Comparing the 
superconducting qubit to the normal conducting qubit 
we notice that here it is e which is flux dependent. 

We can proceed with the discussion of the persistent 
current as in the normal case: In the eigen basis of 
the Hamiltonian the persistent current in the low and 
higher energy states is ±/o = ±dhil/d^ with fin/2 = 
[£;3cos2(7r$^/3>o) + El{l/2 - Ngf]^/^ The operator 
of the persistent current is as for the mesoscopic qubit 
given by 7 = —{2Io/hQ)Hs- The persistent current of 
this superconducting qubit, when it is coupled to a bath 
(external circuit) fluctuates even in the ground state. In 
the experiment of Vion et al. [33] Rabi oscillations are 
reported with an extremely small damping, a Q-factor 
as high as 25'000. This value can be used to estimate 
the coupling constant a and leads to a probability to 
find the system in the excited state of only lO""*. This 
is to small to be of significance when compared to other 
effects, duo the measurement circuit or nearby charge 
traps. The effect we have discussed is most relevant 
in the strong coupling case between system and bath. 
Recent experiments in which the Cooper pair box is 
coupled to a cavity mode [36] offer ways to explore the 
strong coupling limit of interest here. 
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4. Energy fluctuations of the harmonic 
oscillator 



environment causes greater deviation from the Planck 
scale limit, the state loses purity. 



We now consider the entanglement energetics of a 
harmonic oscillator, 



Hs = pV(2m) + (l/2)mtj^g^ 



(20) 



Since there are an infinite number of states, the prob- 
lem is harder. We assume a linear coupling with a har- 
monic oscillator bath. As a consequence the density 
matrix is Gaussian so that environmental information 
is contained in the second moments (g^) and (p^) [5], 



(qIpW) = 



1 



exp 



(21) 

Expectation values of higher powers of Hs are non- 
trivial because q and p do not commute. The purity of 
the density matrix Eq. (21) is 



Trp^ = j dqdq'{q\p\q'){q'\p\q) = 



h/2 



(22) 



The uncertainty relation, \/{cp}{p^ > /i/2, guarantees 
that Trp^ < 1 with the inequality becoming sharp if 
the oscillator is isolated from the environment. As the 




Fig. 5. The probability to measure a harmonic oscillator 
in the ground and first three excited states as a function 
of X and y (see text). The line traces out the behavior 
of the ohmic bath as a function of the coupling in the 
under-damped range. After [3] 



4.1. Energy cumulants 

We calculate the generating function Z{x) = 
{exp{—xHs)) from which we can determine the n**" 
energy cumulant 

((^^r)) = (-)"^ln^(x)|^^^. (23) 
by taking derivatives. Ref. [3] finds 
sinh ex 



2E- 



+ 2 A (cosh £X " 1) 



1 + cosh ex 



(24) 

where e — hcu, 2E = muj^{q'^) + {p^)/m and A = 
{q^) {p^) / . E is the average energy of the oscillator, 
while A > 1 is a measure of satisfaction of the uncer- 
tainty principle. 

Using Eq. (23), the first few harmonic oscillator 
energy cumulants are straightforwardly found via 
Eq. (23), 



((i/f)) = (l/2)[-(eV2) + 4i?' - 2e'A] , 
{{H't)) = -(£/2)[-16£' £^(1 + 12A)] 
{{Hi)) = 48B'' - 4e^S^(l + 12A) 
+ e\(l/%)+2A + QA^]. 



(25) 
(26) 

(27) 



After inserting the mean square values for an ohmic 
bath (see the discussion above Eqs. (29,30)), Eq. (25) 
is identical to Eq. (7). 

4.2. Density matrix 

Alternatively, we now consider the diagonal ma- 
trix elements An analytical expression for the 
density matrix in the energy basis may be found 
by using the wavefunctions of the harmonic oscilla- 
tor, tpn{q) oc ^'^Hn{jq) where 7 = ^ muj/h 
and Hn{x) is the n*'' Hermite polynomial. In the 
energy basis, the density matrix is given by pnm = 
J dqdq'tp^{q){q\p\q')tprn{q')- We first define the dimen- 
sionless variables x — 27^(5^), y = 2{p^)/(7^/i^), and 
D = l + x + y + xy. x and y are related to the major 
and minor axes of an uncertainty ellipse. The isolated 
harmonic oscillator (in it's ground state) obeys two im- 
portant properties: minimum uncertainty (in position 
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and momentum) and equipartition of energy between 

average kinetic and potential energies. The influence 
of the environment causes deviations from these ideal 
behaviors which may be accounted for by introducing 
two new parameters, a = {y — x)/D, b — {xy — 1)/D 
with — 1 < a < 1 and < 6 < 1. The deviation from 
equipartition of energy is measured by o, while the de- 
viation from the ideal uncertainty relation is measured 
by b. We find 



a 



2sn/2 



Pn^b/y/^^] , (28) 



where Pn[z] arc the Lcgcndrc polynomials. The prob- 
ability for the lone oscillator to be measured in an ex- 
cited state clearly decays rapidly with level number. 
These probabilities also reveal environmental informa- 
tion. For example, Pi = 26 and is thus only sensitive 
to the area of the state, while P2 = + 2b^ depends 
on both the uncertainty and energy asj'mmctry. Ad- 
ditionally, if we expand the first density matrix eigen- 
value [5] with respect to small deviations of x and y, we 
recover pn in agreement with our general argument. 

Thus far we have treated x and y as independent vari- 
ables. In reality the environment the system is coupled 
to replaces these variables with two functions of the 
coupling constant. For example, with the ohmic bath 
[5,7] (in the under-damped limit), the variables are 



x{a) 



1/2 a 
1 arctan — , 

,2 V TT ^r^^ 



4a . 



y{a) = (1 - 2a'')x(a) -h — log — , 



(29) 
(30) 



where a is the coupling to the environment in units 

of the oscillator frequency and cUc is a high frequency 
cutoff. This batfi information is shown in Fig. (5) with 
ujc = lOw. The trajectory of the line over the surface 
shows how the probabilities evolve as the coupling a 
is increased from to 1. Other kinds of environments 
would trace out different contours on the probability 
surface. 

4.3. Ground state energy-energy correlations 

The expectation value of observables of the system, 
like the energy Hs , the persistent current I or their mo- 
ments, arc time-independent in the ground state of the 
total system. However, this is not true for two-time cor- 
relations, like the energy-energy correlation function. 



O 



N=200 



N=100 



N=50 



400 



Fig. 6. Energy-energy correlation of a harmonic oscillator 
coupled to a chain of N particles as a function of time. 
The correlation is evaluated in the ground state of the en- 
tire system. The length of the chain is A*" = 50 (bottom), 
N = 100 (middle), and = 200 (top). After Ref. [7]. 



C{t) = -{AHs{t)AHs{0) + AHs{0)AHs{t)). 



(31) 



Here AHs{t) = Ha{t) — {Hs} are the energy fluctua- 
tions away from the average energy. This correlation 
function vanishes if the ground state is a product of a 
system and bath wave function. The correlation is thus 
also a measurement of the degree of entanglement be- 
tween system and bath. 

For the oscillator, Eq. (20), coupled to a linear chain 
of A'^ particles with elongation Xn,n = 1, .., N coupled 
with an energy {l/2)rnhCJh{^n-i—Xn)^ this correlation 
function was calculated by K. E. Nagaev and one of the 
authors [7]. An infinitely long chain generates friction 
proportional to 77 = (mh/m)ujh giving rise to a system 
bath coupling constant a = (m/i/m)(a;(,/a'). The cal- 
culation proceeds by first searching the normal modes 
of the classical problem. The corresponding classical 
problem is then quantized. The energy of the subsys- 
tem and in particular the energy-energy correlation is 
written in terms of the normal modes of the entire sys- 
tem. The correlation is shown in Fig. (6) as a function of 
time for three different chain lengths N = 50, 100, 200 
with Lo/uJh = 1 and nih/m = 0.1. The numbers on the 
vertical axis indicate the initial value of the correlation. 
There is a very rapid initial decay of the correlation 
followed by a damped oscillatory behavior. Since the 
chain is of finite length a partial revival is seen after a 
time it takes a perturbation to travel down the chain 
and back to the oscillator. We emphasize that the re- 
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vival is not complete. The spectrum of the chain alone 

would consists of commensurate frequencies, however, 
due to the presence of the harmonic oscillator the spec- 
trum of the entire system is not commensurate. 

Unlike in the case of the two-state systems we can 
not easily develop a model for the harmonic oscillator 
which connects it to persistent currents. However, the 
conductance in an Aharonov-Bohm geometry in which 
electrons traverse an oscillatory potential has recently 
been discussed by Ratchov et al. [37]. 



5. Discussion 

The energy of the subsystem is an observable which 
illustrates best the distinction between separable and 
entangled ground states. We have shown that pro- 
jective measurements of the system Hamiltonian at 
zero temperature can find the system in higher energy 
states. This is the case if the many-body quantum 
mechanical ground state of system and environment 
are entangled. 

Entanglement assures that the system "knows" 
about its environment and similarly there is system 
information in the environment. If the environment is 
represented as a linear chain of particles the system 
bath interaction can be viewed as a scattering prob- 
lem. System information is then present in the bath 
in the form of the phase of the reflected part of the 
scattering state [38] . 

We emphasize that the system bath interaction is 
not just a question of renormalization. For instance in 
the two state problem discussed here the tunnel matrix 
element A is "renormalized" to Ae//. However, clearly 
the physics we have discussed here is not captured by 
simply replacing in the two-state Hamiltonian the tun- 
nel matrix element by its renormalized value. Both a 
two-state system or a harmonic oscillator (with renor- 
malized mass and frequency) would exhibit a sharp en- 
ergy in their ground state. 

One natural question is how this work is connected 
to the prcscucc or absence of dephasing at zero tem- 
perature. Historically, dephasing has mainly dealt with 
the randomizing of a quantum mechanical phase via 
interaction with some fluctuating variable, such as a 
reduction of the Aharanov-Bohm interference pattern 
from voltage fluctuations, which usually freeze out at 



low temperature. Typically, the off-diagonal density 
matrix elements decay in time, while the diagonal ele- 
ments stay constant. 

A more modern view is to call decoherence any mech- 
anism where one starts with a pure system state, and 
ends with a mixed system state. The ultimate cause 
of this process is simply entanglement of the system 
under observation with other degrees of freedom that 
are unmonitored. Thus, although the entire quantum 
system may be in a pure state, the fact that local mea- 
surements on the subsystem extract only part of the 
information, results in mixed behavior. 

In this sense, there is trivially decoherence at zero 
temperature, unless cither the coupling constant van- 
ishes (so the ground state is separable), or it is possible 
to make measurements on every coupled quantum de- 
gree of freedom, so the purity of the many-body ground 
state is accessible. 

Relaxation into equilibrium is probably the simplest 
possible preparation method of an entangled state. For 
this reason ground state entanglement energetics will 
likely be an important direction of future research. 

This work was supported by the Swiss NSF and by 
the network MaNEP. 
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